import pylab as plt
%matplotlib inline
def read_tab_only(inputfile):
f=open(inputfile,'r').read().split('\n')
time=[]
pitch=[]
for i in range(1,len(f)):
line=f[i]
if line!="":
split=line.split('\t')
#print split
time.append(float(split[0]))
pitch.append(float(split[1]))
#plt.plot(pitch,'gx')
return time,pitch
#'pitch/pitcCHJ000001.tab','procd_pitch_newtrim/CHJ000001_proc.tab'
def plot_one_file(pitch_file,procd_file):
time,pitch=read_tab_only(pitch_file)
ntime,npitch=read_tab_only(procd_file)
plt.figure()
plt.plot(time,pitch,'x')
plt.plot(ntime,npitch,'ro')
plt.title(procd_file)
import random
from os import listdir
pitch_path = 'pitch'
procd_path = 'procd_pitch_newtrim'
onlyfiles = [ f for f in listdir(procd_path) if f.endswith(".tab")]
#print onlyfiles
SEED=249
random.seed(SEED)
sampled_files=random.sample(onlyfiles,50)
print sampled_files
for file_name in sampled_files:
core_name = file_name.split('_')[0]
pitch_file = pitch_path + "/" + 'pitc' + core_name + '.tab'
procd_file = procd_path + '/' + file_name
plot_one_file(pitch_file, procd_file)
import random
from os import listdir
pitch_path = 'pitch'
procd_path = 'procd_pitch_newtrim_sec'
onlyfiles = [ f for f in listdir(procd_path) if f.endswith(".tab")]
#print onlyfiles
SEED=249
random.seed(SEED)
sampled_files=random.sample(onlyfiles,50)
print sampled_files
for file_name in sampled_files:
core_name = file_name.split('_')[0]
pitch_file = pitch_path + "/" + 'pitc' + core_name + '.tab'
procd_file = procd_path + '/' + file_name
plot_one_file(pitch_file, procd_file)
import random
from os import listdir
pitch_path = 'pitch'
procd_path = 'procd_pitch_newtrim_sec_2'
onlyfiles = [ f for f in listdir(procd_path) if f.endswith(".tab")]
#print onlyfiles
SEED=249
random.seed(SEED)
sampled_files=random.sample(onlyfiles,50)
print sampled_files
for file_name in sampled_files:
core_name = file_name.split('_')[0]
pitch_file = pitch_path + "/" + 'pitc' + core_name + '.tab'
procd_file = procd_path + '/' + file_name
plot_one_file(pitch_file, procd_file)
import random
from os import listdir
pitch_path = 'pitch'
procd_path = 'procd_pitch_newtrim_sec_3'
onlyfiles = [ f for f in listdir(procd_path) if f.endswith(".tab")]
#print onlyfiles
SEED=249
random.seed(SEED)
sampled_files=random.sample(onlyfiles,50)
print sampled_files
for file_name in sampled_files:
core_name = file_name.split('_')[0]
pitch_file = pitch_path + "/" + 'pitc' + core_name + '.tab'
procd_file = procd_path + '/' + file_name
plot_one_file(pitch_file, procd_file)
plot_one_file('pitch/pitcXIH000164.tab','procd_pitch_newtrim/XIH000164_proc.tab')